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Abstract 

We study the classical scattering problem of a plane electromagnetic wave incident on the sur- 
face of semi-infinite periodic stratified media incorporating anisotropic dielectric layers with special 
oblique orientation of the anisotropy axes. We demonstrate that an obliquely incident light, upon 
entering the periodic slab, gets converted into an abnormal grazing mode with huge amplitude and 
zero normal component of the group velocity. This mode cannot be represented as a superposition 
of extended and evanescent contributions. Instead, it is related to a general (non-Bloch) Floquet 
eigenmode with the amplitude diverging linearly with the distance from the slab boundary. Re- 
markably, the slab reflectivity in such a situation can be very low, which means an almost 100% 
conversion of the incident light into the axially frozen mode with the electromagnetic energy den- 
sity exceeding that of the incident wave by several orders of magnitude. The effect can be realized 
at any desirable frequency, including optical and UV frequency range. The only essential physical 
requirement is the presence of dielectric layers with proper oblique orientation of the anisotropy 
axes. Some practical aspects of this phenomenon are considered. 
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I. INTRODUCTION 



Electromagnetic properties of periodic stratified media have been a subject of extensive 
research for decades (see, for example, and references therein). Of particular interest 

has been the case of periodic stacks (ID photonic crystals) made up of lossless dielectric com- 
ponents with different refractive indices. Photonic crystals with one-dimensional periodicity 
had been widely used in optics long before the term "photonic crystals" was invented. 

Let us look at the classical problem of a plane electromagnetic wave incident on the 
surface of semi-infinite plane-parallel periodic array, as shown in Fig. ^ The well known 
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FIG. 1: The scattering problem for a semi- infinite periodic layered medium. Si, Sr and St 
are the energy density fluxes of the incident, reflected and transmitted waves, respectively. The 
transmitted wave ^t is a superposition of two Bloch eigenmodes, each of which can be either 
extended or evanescent. Only extended modes can transfer the energy in the z direction. 



effects of the slab periodicity are: (i) the possibility of omnidirectional reflectance when 
the incident radiation is reflected by the slab, regardless of the angle of incidence; (ii) the 
possibility of negative refraction , when the tangential component of the energy flux St of 
transmitted wave is antiparallel to that of the incident wave; (iii) dramatic slowdown of the 
transmitted wave near photonic band edge frequency, where the normal component of the 
transmitted wave group velocity u vanishes along with the respective energy flux St. The 



extensive discussion on the subject and numerous references can be found in 
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All the above effects can occur even in the simplest case of a semi- 



infinite periodic array of two isotropic dielectric materials with different refractive indices, 
for example, glass and air. The majority of known photonic crystals fall into this category. 
The introduction of dielectric anisotropy, however, can bring qualitatively new features to 
electromagnetic properties of periodic stratified media and o pen up new opportunities for 
practical applications (see, for example, a recent publication 14]). One of such phenomena 
is the subject of this work. 

A. The Axially Frozen Mode (AFM) 

Consider a semi-infinite periodic stack with at least one of the constituents being an 
anisotropic dielectric material with oblique orientation of anisotropic axis. A simple example 
of such an array is presented in Fig. El We will show that under certain physical conditions, a 
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FIG. 2: Periodic layered structure with two layers A and B in a primitive cell L. The A layers 
(hatched) are anisotropic with one of the principle axes of the dielectric permittivity tensor e 
making an oblique angle with the normal z to the layers (e xz 7^ 0). The B layers are isotropic. 
The x — z plane coincides with the mirror plane m y of the stack. 

monochromatic plane wave incident on the semi-infinite slab is converted into an abnormal 
electromagnetic mode with huge amplitude and nearly tangential energy density flux, as 
illustrated in Fig. El Such a wave will be referred to as the Axially Frozen Mode (AFM). 
The use of this term is justified because the normal (axial) component u z of the respective 
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FIG. 3: An incident plane wave with unity energy density flux and certain angle of incidence is 
converted into the AFM with huge amplitude, tangential group velocity, and nearly tangential 
energy flux St- The normal components ^ij and (j$T*J of the incident and transmitted waves 
energy flux are comparable in magnitude. 

group velocity becomes vanishingly small, while the amplitude of the AFM can exceed the 
amplitude of the incident plane wave by several orders of magnitude. 

The group velocity u of the AFM is parallel to the semi-infinite slab boundary and, 
therefore, the magnitude of the tangential component (&tJ of the respective energy density 
flux St is overwhelmingly larger than the magnitude of the normal component (j$T^ ■ But, 
although (j^Tj *C (j^Tj j the normal component (j^rj of the energy density flux inside 
the slab is still comparable with that of the incident plane wave in vacuum. This property 
persists even if the normal component u z of the wave group velocity inside the slab vanishes, 
i.e., 

(3r) >0, if «* = (). (1) 

The qualitative explanation for this is that the infmitesimally small value of u z is offset 
by huge magnitude of the energy density W in the AFM. As the result, the product u z W, 
which determines the normal component ($tj of the energy flux, remains finite. The 
above behavior is totally different from what happens in the vicinity of a photonic band 
edge, where the normal component u z of the wave group velocity vanishes too. Indeed, let 
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us introduce the transmittance (r) and the reflectance (p) of a lossless semi-infinite slab 

St) (Sr. 
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In line with Eq. ([Tjl. in the AFM regime the transmittance r remains significant and can be 
even close to unity, as shown in an example in Fig. Ufa). in other words, the incident plane 
wave enters the slab with little reflectance, where it turns into an abnormal AFM with 
infinitesimally small normal component of the group velocity, huge amplitude, and huge 
tangential component of the energy density flux. By contrast, in the vicinity of a photonic 
band edge (at frequencies near u = u>b in Fig. Ufa)), the transmittance of semi- infinite slab 
always vanishes, along with the normal component u z of the wave group velocity. 
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FIG. 4: (a) The transmittance r of periodic semi-infinite slab vs. frequency at fixed direction n of 
the incidence. At the frequency ojq of the AFM, r is close to unity, which implies that the incident 
wave almost completely gets converted into the AFM. (b) The respective axial dispersion relation 
u){k z ) at fixed (n x ,n y ) from Eq. (12). At k z = ko and to = u>q this spectral branch develops a 
stationary inflection point (16) associated with the AFM regime. cj& is the edge of the frequency 
band for a given (n x ,n y ). The values of u and k are expressed in units of c/L and 1/L, respectively. 

It turns out that at a given frequency ujq the AFM regime can occur only for a special 
direction n of the incident plane wave propagation 



n = n (u ) 



(3) 
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This special direction of incidence always makes an oblique angle with the normal z to the 
layers. To find no for a given ujq or, conversely, to find o>o for a given fto, one has to solve 
the Maxwell equations in the periodic stratified medium. This problem will be addressed in 
Section 3. In Section 2 we consider the relation between the AFM regime and the singularity 
of the electromagnetic dispersion relation responsible for such a peculiar behavior. If the 
frequency uo and the direction of incidence n do not match explicitly as prescribed by Eq. 

(0) , the AFM regime will be somewhat smeared. 

B. The vicinity of the AFM regime 

Let {z) be the transmitted electromagnetic field inside the semi-infinite slab (the 
explicit definition of \&t (z) is given in Eqs. (}3*H|) and (jHHJ)). It turns out that in the vicinity of 
the AFM regime, (z) is a superposition of the extended and evanescent Bloch eigenmodes 

^ T (z) = V ex (z) + m ev (z), z> 0, (4) 

where ^ ex (z) is an extended mode with u z > 0, and \E% (z) is an evanescent mode with 
Qk z > 0. As shown in an example in Fig. both the contributions to ( z ) have huge and 
nearly equal and opposite values near the slab boundary, so that their superposition (jU) at 
z = is small enough to satisfy the boundary condition (|90|). As the distance z from the 
slab boundary increases, the evanescent component *f? ev (z) decays exponentially, while the 
amplitude of the extended component ^ ex (z) remains constant and huge. As the result, the 
field amplitude (z)\ 2 reaches its huge saturation value \^ e x\ 2 at a certain distance from 
the slab boundary (see Eqs. O, (flip and AMD). 

When the direction of incidence n tends to its critical value no for a given frequency uoq, the 
respective saturation value |\l/ex| 2 of the AFM amplitude (z)\ 2 diverges as \n — no\~ 2 ^ 3 . 
Conversely, when the frequency uo tends to its critical value ujq for a given direction of 
incidence no, the saturation value of the AFM amplitude diverges as \u — ujq\~ 2 ^ . In the 
real situation, of course, the AFM amplitude will be limited by such physical factors as: 

(1) nonlinear effects, (ii) electromagnetic losses, (iii) structural imperfections of the periodic 
array, (iv) finiteness of the slab dimensions, (v) deviation of the incident radiation from a 
perfect plane monochromatic wave. 

Fig. El gives a good qualitative picture of what really happens in the vicinity of the AFM 
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a) Resulting field b) Extended component c) Evanescent component 
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FIG. 5: Destructive interference of the extended and evanescent components of the resulting electro- 
magnetic field (4) inside semi-infinite slab in close proximity of the AFM regime: (a) the amplitude 
( z )\ 2 °f the resulting field, (b) the amplitude \^ ex (z)\ 2 of the extended contribution, (c) the 
amplitude |\l/ e ,; (z)| 2 of the evanescent contribution. The amplitude of the incident wave is 

unity. The distance z from the slab boundary is expressed in units of L. 

regime. Consider a wide monochromatic beam of frequency uo incident on the surface of semi- 
infinite photonic slab. The direction of incidence n \\ Si is chosen so that the condition (jHJ) 
of the AFM regime is satisfied at u = ujq- As frequency us tends to u>o from either direction, 
the normal component u z of the transmitted wave group velocity approaches zero, while the 
tangential component u± remains finite 

u z ~ \oj — co>o| 2 ^ 3 — ► 0, u± — > uq as io — - > c^o- (5) 

This relation together with the equality 

7T U z 

Vt — arctan — (6) 

2 u± 

involving the refraction angle 9t, yield 

— — 6t ~ \u) — uj \ 2 ^ 3 ^0 as u-> lu . (7) 

Hence, in the vicinity of the AFM regime, the transmitted (refracted) electromagnetic wave 
can be viewed as a grazing mode. The most important and unique feature of this grazing 
mode directly relates to the fact that the transmittance r of the semi-infinite slab remains 
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FIG. 6: Incident and transmitted (refracted) waves in the vicinity of the AFM regime. The reflected 
wave is not shown. 9i and 9t are the incidence and refraction angles, Si and St are the energy 
density fluxes of the incident and transmitted waves. Both the energy density and the energy 
density flux in the transmitted wave are much larger than the respective values in the incident 
wave. However, the total power transmitted by the refracted wave is smaller by factor r, due to 
much smaller cross-section area of the nearly grazing transmitted wave. 

finite even at uj = u (see, for example, Fig. Ufa)). Indeed, let Aj and A T be the cross-section 
areas of the incident and transmitted (refracted) beams, respectively. Obliviously, 

At _ cos6> T 

A! ~ COS0/ 1 ' 

Let us also introduce the quantities 

Uj = AjSj, U T = A T S T} (9) 



S 



where Sj and St are the energy density fluxes of the incident and transmitted waves. Uj 
and Ut are the total power transmitted by the incident and transmitted (refracted) beams, 
respectively. The expressions (JBJ) and (JOJ) imply that 

U T S t cos 9 t (S t ) z 
Ui Si cos 1 {bi) z 

which is nothing more than a manifestation of the energy conservation law. Finally, Eq. 
([EH) , together with the formula (JZj), yield 

C C C ° S ^ I 1-2/3 /„n 

St = rbi — ~ a; — a; — > oo as c<j — > u; - (H) 

cos # r 

where we have taken into account that tSj cos 8j is limited (of the order of magnitude of 
unity) as to — > ujq. By contrast, in the vicinity of the photonic band edge the transmittance 
t of the semi-infinite slab vanishes along with the energy density flux St of the transmitted 
(refracted) wave. 

The expressions (JZJ) and show that in the vicinity of the AFM regime, the transmitted 
wave behaves like a grazing mode with huge and nearly tangential energy density flux St 
and very small (compared to that of the incident beam) cross-section area At, so that the 
total power Ut = AtSt associated with the transmitted wave cannot exceed the total power 
Ui of the incident wave: Ut = tUj < Uj. 

The above qualitative consideration is only valid on the scales exceeding the size L of 
the unit cell (which is of the order of magnitude of c/u) and more importantly, exceeding 
the transitional distance I = (S/Cg^) -1 from the slab boundary where the evanescent mode 
contribution to the resulting electromagnetic field \I/t (z) is still significant. The latter 
means that the width of both the incident and the refracted beams must be much larger 
than I. If the above condition is not met, we cannot treat the transmitted wave as a beam, 
and the expressions (J7J) through (fTTj) do not apply. Instead, we would have to use the 
explicit electrodynamic expressions for ^ (z), such as the asymptotic formula (jlUlj) . Note 
that if the direction n of the incident wave propagation and the frequency u exactly match 
the condition Q for the AFM regime, the transmitted wave *&t ( z ) does not reduce to 
a superposition (JH) of canonical Bloch eigenmodes. Instead, the AFM is described by a 
general Floquet eigenmode \l/oi {%) from Eq. (JKUj) . which diverges inside the slab as z, until 
the nonlinear effects or other limiting factors come into play. The related mathematical 
analysis is provided in Sections 3 and 4. 
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In some respects, the remarkable behavior of the AFM, is similar to that of the frozen 
mode related to the phenomenon of electromagnetic unidirectionality in nonreciprocal mag- 



netic photonic crystals 
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16(. In a unidirectional photonic crystal, electromagnetic radia- 



tion of a certain frequency ujq can propagate with finite group velocity u\\z only in one of the 
two opposite directions, say, from right to left. The problem with the electromagnetic uni- 
directionality, though, is that it essentially requires the presence of magnetic materials with 
strong circular birefringence (Faraday rotation) and low losses at the frequency range of in- 
terest. Such materials are readily available at the microwave frequencies, but at the infrared 
and optical frequency ranges, finding appropriate magnetic materials is highly problematic. 
Thus, at frequencies above 10 12 Hz, the electromagnetic unidirectionality along with the 
respective nonreciprocal magnetic mechanism of the frozen mode formation may prove to be 
impractical. By contrast, the occurrence of AFM does not require the presence of magnetic 
or any other essentially dispersive components in the periodic stack. Therefore, the AFM 
regime can be realized at any frequencies, including the infrared, optical, and even ultraviolet 
frequency ranges. The only essential physical requirement is the presence of anisotropic di- 
electric layers with proper orientation of the anisotropy axes. An example of such an array 
is shown in Fig. El 

In Section 2 we establish the relation between the phenomenon of AFM and the electro- 
magnetic dispersion relation of the periodic layered medium. This allows us to formulate 
strict and simple symmetry conditions for such a phenomenon to occur, as well as to find 
out what kind of periodic stratified media can exhibit the effect. Relevant theoretical anal- 
ysis based on the Maxwell equations in stratified media is carried out in Sections 3 and 4. 
Finally, in Section 5 we discuss some practical aspects of the phenomenon. 



II. DISPERSION RELATION WITH THE AFM 



Now we establish the connection between the phenomenon of AFM and the electromag- 
netic dispersion relation u (jtj , k = (k x , k y , k z ) of the periodic stratified medium. In a 
plane-parallel stratified slab, the tangential components (k x , k y ) of the Bloch wave vector k 
always coincide with those of the incident plane wave in Figs. El and El while the normal 
component k z is different from that of the incident wave. To avoid confusion, in further 
consideration, the z component of the Bloch wave vector k inside the periodic slab will be 
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denoted as k without the subscript z, namely 



Inside periodic stack: k = (k x , k y , k). 

The value of k is found by solving the Maxwell equations in the periodic stratified medium 
for given u and (k x , k y ); k is defined up to a multiple of 2-n / L, where L is the period of the 
layered structure. 

Consider now the frequency u as function of k for fixed (k x ,k y ). A typical example of 
such a dependence is shown in Fig. |7[a). A large gap at the lowest frequencies is determined 
by the value of the fixed tangential components (k x , k y ) of the quasimomentum k. This gap 
vanishes in the case of normal incidence, when k x = k y = 0. An alternative and more 
convenient representation for the dispersion relation is presented in Fig. E^b), where the 
plot of u(k) is obtained for fixed (n x ,n y ) based on 



(n x ,n y ) = (ck x /u,ck y /u). 



(12) 



The pair of values (n x ,n y ) coincide with the tangential components of the unit vector n 
defining the direction of the incident plane wave propagation. The dependence u(k) for 
fixed (n x ,n y ) or for fixed (k x , k y ) will be referred to as the axial dispersion relation. 

Suppose that for k = k and uj = too = oj (koj j one °f the spectral branches u (k) develops 
a stationary inflection point for given (k x , k y ) = (k 0x , k 0y ), i.e., 
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(14) 



in Eq. ()13j) is the axial component of the group velocity, which vanishes at k — ko- Observe 
that 

/ ^'^ X f 9uj\ 

and Ui, = I ■77; — I , (15) 



u , 



dto 

dk x j f , 



dk. 



y ' k.kx 



representing the tangential components of the group velocity, may not be zeros at k = k$. 

Notice that instead of (|13)1. one can use another definition of the stationary inflection 
point 
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FIG. 7: The axial dispersion relation of anisotropic periodic stack in Fig. 2: (a) u(k z ) for 
fixed values (k x , ky) of the tangential components of quasimomentum k; (b) ui{k z ) for fixed values 
(n x ,n y ), defining the direction of incidence. In the case of normal incidence, there would be no 
difference between (a) and (b). 

The partial derivatives in Eqs. (|T6|) are taken at constant (n x ,n y ), rather than at constant 
(k x , k y ). Observe that the definitions (|T3j) and (fTo*}) are equivalent, and we will use both of 
them. 

In Fig. IH^b) we reproduced an enlarged fragment of the upper spectral branch of the axial 
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dispersion relation in Fig. W[b). For the chosen (n x ,n y ), this branch develops a stationary 
inflection point (fTSJl at uj = ujq and k = ko- The extended Bloch eigenmode with uj = ujq 
and k = k , associated with the stationary inflection point, turns out to be directly related 
to the axially frozen mode ( AFM) . 

In Sections 3 and 4, based on the Maxwell equations, we prove that the singularity 
(or, equivalently, (|T3"j) ) indeed leads to the very distinct AFM regime in the semi- infinite 
periodic stack. We also show that a necessary condition for such a singularity and, therefore, 
a necessary condition for the AFM existence is the following property of the axial dispersion 
relation of the periodic stack 

uj (k x , k y , k) ^ uj (k x , k y , —k) or, equivalently, uj (n x , n y , k) ^ to {n x , n y , —k) (17) 

This property will be referred to as the axial spectral asymmetry. Evidently, the axial 
dispersion relations presented in Fig. satisfy this criterion. Leaving the proof of the 
above statements to Section 3, let us look at the constraints imposed by the criterion (fTTjl 
on the geometry and composition of the periodic stack. 

A. Conditions for the axial spectral asymmetry 

First of all, notice that a periodic array would definitely have an axially symmetric dis- 
persion relation 

uj (k x , k y , k) = uj {k x , k y , —k) or, equivalently, uj (n x , n y , k) = uj {n x , n y , —k) (18) 

if the symmetry group G of the periodic stratified medium includes any of the following two 
symmetry operations 

m z , 2' z = 2 Z x R, (19) 

where m z is the mirror plane parallel to the layers, 2 Z is the 2-fold rotation about the z 
axis, and R is the time reversal operation. Indeed, since 2 Z 
R {k x , k y , k) = (—k x , —k y , —k), we have 

2 2 (k x , ky, k) = {k x , ky, —k) , 

which implies the relation (|18jl for arbitrary (k x , k y ). The same is true for the mirror plane 

m z 

m z {k x , ky, k) = {k x , ky, k) . 
13 



Consequently, a necessary condition for the axial spectral asymmetry (fTTj) of a periodic stack 
is the absence of the symmetry operations ([19)1 . i.e., 

m z i G and 2' z <£ G. (20) 

In reciprocal (nonmagnetic) media, where by definition, R e G, instead of Eq. ("2*0*)) one can 
use the following requirement 

m z (£G and 2 Z <£ G. (21) 
Note, that the axial spectral symmetry ("I~8j) is different from the bulk spectral symmetry 

u) (k x , ky, k) = lu ( k X} ky } k) (22) 

For example, the space inversion / and/or the time reversal R, if present in G, ensure the 
bulk spectral symmetry (|22jl . but neither / nor R ensures the axial spectral symmetry (|18jl. 

1. Application of the criterion h21\) to deferent periodic stacks. 

The condition ("2~T~) for the axial spectral asymmetry imposes certain restrictions on the 
geometry and composition of the periodic stratified medium, as well as on the direction of 
the incident wave propagation. 

a. Restrictions on the geometry and composition of periodic stack. First of all, observe 
that a common periodic stack made up of isotropic dielectric components with different 
refractive indices, always has axially symmetric dispersion relation (fTHj) . no matter how 
complicated the periodic array is or how many different isotropic materials are involved. To 
prove this, it suffices to note that such a stack always supports the symmetry operation 2 Z . 

In fact, the symmetry operation 2 Z holds in the more general case when all the layers are 
either isotropic, or have a purely in-plane anisotropy 
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(23) 



Obviously, the in-plane anisotropy ([2*3*1) does not remove the symmetry operation 2 Z and, 
therefore, the property (|18j) of the axial spectral symmetry holds in this case. Thus, we 
can state that in order to display the axial spectral asymmetry, the periodic stack must 
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include at least one anisotropic component, either uniaxial or biaxial. In addition, one of 
the principle axes of the respective dielectric permittivity tensor e must make an oblique 
angle with the normal to the layers, which means that at least one of the two components 
e xz and e yz of the respective dielectric tensor must be nonzero. 

The above requirement gives us a simple and useful idea on what kind of periodic stratified 
media can support the axial spectral asymmetry and the AFM regime. But this is not a 
substitute for the stronger symmetry criterion (|2"UJl or (|21|). For example, although the 
periodic stack in Fig. |S] includes the A layers identical to those in Fig. El this stack does 
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FIG. 8: Periodic stack composed of anisotropic layers A\ and A2, which are the mirror images of 
each other, and isotropic layers B. This stack has axially symmetric dispersion relation and does 
not support the AFM regime. This is true even if the B layers are removed. 

not meet the criterion ()2U|) for the axial spectral asymmetry. Indeed, the stack in Fig. |H] 
supports the mirror plane m z , which, according to the expression ()19p . ensures the axial 
spectral symmetry. 

b. Restriction on the direction of incident wave propagation Consider now an impor- 
tant particular case k x — k y — of the normal incidence. The criterion (|T7|) reduces now to 
the simple requirement 

u (k) ^ u) (—k) , where k — (0, 0, k) , (24) 

of the bulk spectral asymmetry, which is prohibited in nonmagnetic photonic crystals due 
to the time reversal symmetry. Therefore, in the nonmagnetic case, we have the following 
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additional condition for the axial spectral asymmetry 

k± = yjkl + kl + 0, (25) 

implying that the AFM cannot be excited in a nonmagnetic semi-infinite stack by a normally 
incident plane wave, i.e., the incident angle must be oblique. 

Conditions (}2~T]) and ()25l) may not be necessary in the case of nonreciprocal magnetic 
stacks (see the details in [15]). But as we mentioned earlier, at frequencies above 10 12 
Hz, the nonreciprocal effects in common nonconducting materials are negligible. Therefore, 
in order to have a robust AFM regime in the infrared or optical frequency range, we must 
satisfy both requirements (|21|) and (|25|). regardless of whether or not nonreciprocal magnetic 
materials are involved. 

As soon as the above conditions are met, one can always achieve the AFM regime at 
any desirable frequency uj within certain frequency range Aou. The frequency range Auj is 
determined by the stack geometry and the dielectric materials used, while a specific value 
of uj within the range can be selected by the direction n of the light incidence. 



B. Periodic stack with two layers in unit cell 

The simplest and the most practical example of a periodic stack supporting the axial 
spectral asymmetry (fT7|) and, thereby, the AFM regime, is shown in Fig. |21 It is made up of 
anisotropic A layers alternated with isotropic B layers. The respective dielectric permittivity 
tensors are 
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For simplicity, we assume 

Va = I l b = I- ( 27 ) 
The stack in Fig. Elhas the monoclinic symmetry 

2 y /m y (28) 

with the mirror plane m y normal to the y - axis. Such a symmetry is compatible with the 
necessary condition ()21|) for the AFM existence. But as we will see below, the symmetry 
f!28|) imposes additional constraints on the direction n of the incident wave propagation. 
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In Fig. El we show the axial dispersion relation ui (k) of this periodic array, computed for 
four different directions (n x , n y ) of incident wave propagation. These four cases cover all the 
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FIG. 9: Axial dispersion relation u) (k) for fixed (n x ,n y ) for the periodic array in Fig. 2. The AFM 
regime can occur only if n x ^ and n y / (the case (d)). 



possibilities, different in terms of symmetry. 

In the case (a) of normal incidence, when n x = n y = 0, the dispersion relation is axially 
symmetric, as must be the case with any reciprocal periodic stratified medium (see the 
explanation after Eq. 
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In the case (b), when n x = and n y ^ 0, the two necessary conditions (|21|) and ()25|) 
for the axial spectral asymmetry are met. Yet, those conditions prove not to be sufficient. 
Indeed, if n x = 0, either of the symmetry operations 

2 y and m' y = m y x R (29) 

imposes the relation 

u(0,k y ,k) =u{0,k y ,-k) (30) 

which implies the axial spectral symmetry. Neither stationary inflection point, nor AFM 
can occur in this case. 

In the case (c), when n x ^ and n y = 0, the situation is more complicated. The 
quasimomentum k lies now in the x — z plane, which coincides with the mirror plane m y . 
Therefore, every Bloch eigenmode (z) can be classified as a pure TE or pure TM mode, 
depending on the (z) parity with respect to the mirror reflection m y 

for TE mode m y ^^ (z) = (z) ; for TM mode: ra^g (z) = ^ (z) . (31) 

The TE modes have axially symmetric dispersion relation 

for TE modes: u (k x , 0, k) = lu (k x , 0, —k) . (32) 

Indeed, the component e xz of the dielectric tensor ea does not affect the TE modes, because 
in this case the electric component E(r, t) of the electromagnetic field is parallel to the y 
axis. As a consequence, the axial dispersion relation of the TE spectral branches is similar 
to that of the isotropic case with e xz = 0, where it is always symmetric. By contrast, for the 
TM modes we have E(r, t) _L y. Therefore, the TM modes are affected by e xz and display 
axially asymmetric dispersion relation 

for TM modes: u (k x , 0, k) ^ uj (k x , 0, — k) , (33) 

as seen in Fig. El^c). We wish to remark, though, that the equality (J32|) cannot be derived 
from symmetry arguments only. The axial spectral symmetry of the TE modes is not exact 
and relies on the approximation (|2"Tj) for the magnetic permeability of the A layers. On 
the other hand, the fact that the spectral branches have different parity (}3*Tj) with respect 
to the symmetry operation m y , implies that none of the branches can develop a stationary 
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inflection point (see Eq. and explanations thereafter). Thus, in the case n y = 0, in 

spite of the axial spectral asymmetry, the AFM regime cannot occur either. 

Finally, in the general case (d), when n x 7^ and n y 7^ 0, all the spectral branches display 
the property (|17j) of the axial spectral asymmetry. In addition, the Bloch eigenmodes now 
are of the same symmetry (i. e., belong to the same irreducible representation of the wave 
vector symmetry group) and are neither TE, nor TM. This is exactly the case when the 
AFM regime can be achieved at some frequencies by proper choice of the incident angle. 
For instance, if we impose the equality n x = n y and change the incident angle only, it turns 
out that every single spectral branch at some point develops a stationary inflection point 
and, thereby, displays the AFM at the respective frequency. If we want the AFM at a 
specified frequency u)q, then we will have to adjust both n x and n y . 

III. ELECTRODYNAMICS OF THE AXIALLY FROZEN MODE 
A. Reduced Maxwell equations 

We start with the classical Maxwell equations for time-harmonic fields in nonconducting 
media 

V x E (f) = i-B (f) , V x H (f) = -i-B (f) , (34) 
c c 

where 

D (r) = e (f) E (f) , B (f) = fi (f) H (f) , (35) 

In a lossless dielectric medium, the material tensors e (r) and fi (r) are Hermitian. In a 
stratified medium, the tensors e (r) and fx (r) depend on a single Cartesian coordinate z, 
and the Maxwell equations ()34j) can be recast as 

V x E (f) = i-jx (z) H (f) , V x H (f) = -%-e (z) E (f) . (36) 

c c 

Solutions for Eq. (|36j) are sought in the following form 

E (f) = e ^+ k vy)E (z) , H (f) = e i{ - k * x+k ^H (z) . (37) 
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The substitution (|37jl transforms the system of six linear equation into a system of four 
linear differential equations 



d z V(z) = i-M{z)V{z), V(z) 
c 



The explicit expression for the Maxwell operator M (z) is 



£ y (z) 
H x (z) 

Hy(z) 



(38) 
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M 2 i M 22 
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The Cartesian components of the material tensors e and fx are functions of z and (in dis- 
persive media) w. The reduced Maxwell equation (J3H|) should be complemented with the 
following expressions for the z components of the fields 



(40) 



E z = (—n x H y + n y H x — e* lz E x — e* 23 E y ) e zz l 
H z = (n x E y - n y E x - li* 13 H x - Lt* 23 H y ) 

where ( defined in Eq. ffTSl). 

Notice that in the case of normal incidence, the Maxwell operator is drastically simplified 

Mn = M 22 = 0, for n x = n y = 0. (41) 

This is the case we dealt with in when considering the phenomenon of electromagnetic 
unidirectionality in nonreciprocal magnetic photonic crystals. By contrast, the objective of 
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this Section is to show how the terms Mu and M22, occurring only in the case of oblique 
incidence, can lead to the phenomenon of AFM, regardless of whether or not the nonrecip- 
rocal effects are present. Note that M%% and M22 are also nonzero in materials with linear 
magnetoelectric effect (see, for example, Ref. 15 and references therein), but we are not 
considering here such an exotic situation. 

Importantly, the 4x4 matrix M (z) in Eq. (|39|) has the property of J - Hermitivity 
defined as 

(JM) t = JM (42) 

where 



j = .r 1 



0001 
00-10 
0-100 
1000 



(43) 



Different versions of the reduced Maxwell equation ()38j) can be found in the extensive 



literature on electrodynamics of stratified media (see, for example, 17j,ll8|, 3], and references 
therein). For more detailed studies of J - Hermitian and J - unitary operators see 



B. The transfer matrix 



The Cauchy problem 



d z V(z)=i-M(z)*(z), *(z )=% 
c 



(44) 



for the reduced Maxwell equation (j38|) has a unique solution 



*(z) = T{z,Zo)*(zo) 



(45) 



where the 4x4 matrix T (z, Zq) is so-called transfer matrix. From the definition (|43|) of the 
transfer matrix it follows that 



T (z, Zq) = T (z, z') T (zf, z ) , T (z, z ) =T 1 (z , z) , T (z, z) = I. 



(46) 



The matrix T (z, z ) is uniquely defined by the following Cauchy problem 



.U) 



d z T (z, Zq) = i-M (z) T (z, Zq) , T (z, z) = I. 
c 



(47) 
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The equation (|47|) . together with J - Hermitivity (jH)) of the Maxwell operator M (z), imply 
that the matrix T (z, zq) is J - unitarity, i.e., 

T^(z,z ) = JT- l (z,z )J. (48) 

(see the proof in Appendix 1). The J - unitarity ()48|) of the transfer matrix imposes strong 
constraints on its eigenvalues (see Eq. It also implies that 

\detT(z,z )\ = l. (49) 

The transfer matrix Ts of a stack of layers is a sequential product of the transfer matrices 
T m of the constitutive layers 

Ts = Y[T m (50) 

m 

If the individual layers are homogeneous, the corresponding single-layer transfer matrices 
T m are explicitly expressed in terms of the respective Maxwell operators M m 

T m = exp (iz m M m ) (51) 

where z m is the thickness of the m-th layer. The explicit expression for M m is given by 
(j2HJ)- Thus, formula (|5Uj) . together with (J5TJ) and (|3*3|). gives us an explicit expression for 
the transfer matrix T5 of an arbitrary stack of anisotropic dielectric layers. Ts is a function 
of (i) the material tensors e and p, in each layer of the stack, (ii) the layer thicknesses, (hi) 
the frequency u, and (iv) the tangential components (k x ,k y ) = (n x ui / c, n y u> / c) of the wave 
vector. 

Consider the important particular case of normal wave propagation. Using Eq. (|HT|) and 
the explicit expression (J39|) for the Maxwell operator, one can prove that 

det (Ts) = 1, for n x = n y = 0. (52) 

Additional information related to the transfer matrix formalism can be found in 
and references therein. 

C. Periodic arrays. Bloch eigenmodes. 

In a periodic layered structure, all material tensors, along with the J - Hermitian matrix 
M(z) in Eq. (JHEJ), are periodic functions of z 

M(z + L) = M(z) (53) 
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where L is the length of a primitive cell of the periodic stack. By definition, Bloch solutions 
(z) of the reduced Maxwell equation (|S5|l with the periodic operator M(z) satisfy 

(z + L) = e lkL V k (z) (54) 

The definition (}4*5|) of the T - matrix together with Eq. (}5lj) give 

V k (z + L) = T (z + L, z) V k (z) = e tkL V k (z) . (55) 

Introducing the transfer matrix of a primitive cell 

T L = T(L,0) (56) 

we have from Eq. (JH5J) 

T L $ fc = e ifcL $ fc , where $ fc = * fc (0) . (57) 

Thus, the eigenvectors of the transfer matrix Tl of the unit cell are uniquely related to the 
eigenmodes of the reduced Maxwell equation (J58j) through the relations 



$ fcl = *fcx (0) , = ^ fc2 (0) , $ fe = fc3 (0) , $ fc4 = ^ fe4 (0) (58) 
The respective four eigenvalues 

Xi = e ihL } i = 1,2,3,4 (59) 

of Tl are the roots of the characteristic equation 

F (X) = 0, where F (X) = det (t l - xfj = X 4 + P 3 X 3 + P 2 X 2 + P X X + 1. (60) 

For any given u and (k x ,k y ), the characteristic equation defines a set of four values 
{Xl, X 2 , X 3 , X 4 }, or equivalently, {^1,^2,^3,^4}. Real k correspond to propagating Bloch 
waves (extended modes), while complex k correspond to evanescent modes. Evanescent 
modes are relevant near photonic crystal boundaries and other structural irregularities. 

The J-unitarity (J4*H|) of Tl imposes the following restriction on the eigenvalues (J59j) for 
any given u and (k x ,k y ) 

{ki} = {k*}, i = 1,2, 3, 4. (61) 
In view of the relation (jblj). one has to consider three different situation. The first possibility 

k\ = fc*, k<z = k 2 , k 3 = fcg, ki = k\ (62) 
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relates to the case of all four Bloch eigenmodes being extended. The second possibility 

k\ = fc*, k 2 = k 2 , &3 = k\, where k 3 ^ k 3 , k± ^ k* Al (63) 
relates to the case of two extended and two evanescent modes. The last possibility 

k\ = k 2 , k 3 = kl, where k\ ^ k*, k 2 ^ k 2 , k 3 ^ k^, k± ^ k\ (64) 

relates the case of a frequency gap, when all four Bloch eigenmodes are evanescent. 
Observe that the relation 

fci + k 2 + k 3 + k 4 = 



valid in the case of normal incidence (see Refs. 



15 



161]). may not apply now. 



1. Axial spectral symmetry 

Assume that the transfer matrix Tl is similar to its inverse 

T L = XJ- X T L X XJ (65) 

where U is an invertible 4x4 matrix. This assumption together with the property (j48jl of 
J-unitarity, imply the similarity of Tl and Tjf 

T l = V~ x T]y, where V = JU. (66) 

This relation imposes additional restrictions on the eigenvalues (|59|) for a given frequency uo 
and given (k x ,k y ) 

{h} = {-hi}, i = 1,2,3,4. (67) 

The relation (|6T|) is referred to as the axial spectral symmetry, because in terms of the 
corresponding axial dispersion relation, it implies the equality |T8|) for every spectral branch. 

If the sufficient condition ()65)1 for the axial spectral symmetry is not in place, then we 
can have for a given u and (k x , k y ) 

{h}^{~h}, % = 1,2,3,4 (68) 
which implies the axial spectral asymmetry ()17j) . 
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D. Stationary inflection point 



The coefficients of the characteristic polynomial F (X) in Eq. ()60|) are functions of uo 
and (k x , k y ). Let Fq (X) be the characteristic polynomial at the stationary inflection point 
(|16p . where u = Uq and (k x ,k y ) = (A<):r? ^oy)- The stationary inflection point (jTSj) can also 
be defined as follows 

F (X) = 0, Fq (X) = 0, Fq (X) = 0, Fq (X) ^ 0. (69) 

This relation requires the respective value of X = exp {ik^L) to be a triple root of the 
characteristic polynomial F (X) implying 

F (X) = (X - Xi) (X - X ) 3 = 0. (70) 

A small deviation of the frequency u from its critical value uq changes the coefficients of 
the characteristic polynomial and removes the triple degeneracy of the solution X 

1/3 



2m/ 3 



(71) 



Y X ~ fiV3 f 9F o/ duJ \ (., \l/3 t t _ 1 p 2^/3 _-2 

or, in terms of the axial quasimomentum k 

k-h~ 6 1 / 3 {^p) V3 & e = 1, e 2m/3 , e- 2 -/ 3 (72) 



where 



in 



The three solutions ([72)1 can also be rearranged as 



> 0. (73) 

k=kg 



k ex *k Q + &/* K")- 1/3 (^-^o) V3 , 

^ « ^0 + \ (6) 1/3 K')- 1/3 ( W - ^ ) 1/3 + ^6V3 (^)-V3 | W _ Wo |l/3 ? (74) 

W«*b+| (6) 1/3 K')- 1/3 ( W - ^ ) 1/3 - ^6V3 KT V3 | W _ Wo |i/a . 

The real k ex in (|74|) relates to the extended mode \l/ ex (2;), with u z = at uj = ojq. The 
other two solutions, k ev and k EV = k* ev) correspond to a pair of evanescent modes ty ev (z) 
and ^ ev (z) with positive and negative infinitesimally small imaginary parts, respectively. 
Those modes are truly evanescent (i.e., have ^k ^ 0) only if oj ^ oj , but it does not mean 
that at oj — ojq, the eigenmodes ^f ev (z) and *f?EV i z ) become extended. In what follows we 
will take a closer look at this problem. 
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1. Eigenmodes at the frequency of AFM 



Consider the vicinity of stationary inflection point (|T3j). As long as w ^ uq, the four 
eigenvectors (J58|) of the transfer matrix Tl comprise two extended and two evanescent Bloch 
solutions. One of the extended modes (say, corresponds to the non-degenerate real 
root X\ = e tklL of the characteristic equation (|60p. This mode has negative axial group 
velocity u z (ki) < and, therefor, is of no interest for us. The other three eigenvectors 
of Tl correspond to three nearly degenerate roots (|7T|). As uj approaches ojq, these three 
eigenvalues become degenerate, while the respective three eigenvectors $fc 2 ,$fe 3 , and 
become collinear 



ai$ fco , $ fc3 -> a 2 $ ko , -> a3$fc , as w 



(75) 



The latter important feature relates to the fact that at uj = Uq, the matrix Tl has a nontrivial 
Jordan canonical form 



U~ l T L U 



X x 
A 1 
X 1 
X 



at w = ciJo 



(76) 



and, therefore, cannot be diagonalized. It is shown rigorously in , that the very fact that 
the Tl eigenvalues display the singularity (|7T|). implies that at uj — u , the matrix Tl has 
the canonical form ([76)1 . In line with (|75j). the matrix Tl from Eq. (|75j) has only two (not 
four!) eigenvectors: 

1. = \l/fc 1 (0), corresponding to the non-degenerate root X% and relating to the ex- 
tended mode with u z < 0; 

2. $fc = \l/fc (0), corresponding to the triple root X and related to the AFM. 

The other two solutions of the Maxwell equation (J38|) at uj = ujq are general Floquet 
eigenmodes, which do not reduce to the canonical Bloch form (|54p. Yet, they can be related 
to \l/fco ( z )- Indeed, following the standard procedure (see, for example, |2lL l22j). consider 
an extended Bloch solution \I/ fe (z) of the reduced Maxwell equation 



.UJ 



L^ fc (z) = 0, where L = d z - i—M (z) 



(77) 
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where both operators M (z) and L (z) are functions of uo and (k x , k y ). Assume now that the 
axial dispersion relation u (k) has a stationary inflection point ()13j) at k — k$. Differentiating 
Eq. (|77|) with respect to k at constant (k x , k y ) gives, with consideration for Eq. (JT3J), 

W k ^! k (z) = 0, W 2 kk ^! k (z) = 0, at k = k . 

This implies that at k = k , both functions 

v& 01 (z) = d k V k (z)\ k=ko , and * 02 (z) = & kk V k (z)\ k=kQ (78) 

are also eigenmodes of the reduced Maxwell equation at u — ojq. Representing ^ k (z) in the 
form 

# fc ( z ) = fa ( z ) e ikz , where ^ k (z + L) = ^ k (L) , %k = (79) 
and substituting Eq. (fT9*j) into (fTSj) we get 

*oi (z) = *fc (*) + ^*fco (*) > ( 80 ) 
^02 (z) = % (z) + iz$ k0 (z) - z 2 m h) (z) , (81) 



' k=k(j ^ 



ikgz 



where 

* fco (z) = (d k fa (z)) k=ko e *»* and % o (z) = (d 2 kk ^ k (z)), 

are auxiliary Bloch functions (not eigenmodes). 

To summarize, at the frequency c^o of AFM, there are four solutions for the reduced 
Maxwell equation ([55)1 

* kl (z) , ^ fc0 (z) , * i (z) , ^02 (*) (82) 

The first two solutions from (|H2*|) are extended Bloch eigenmodes with u z < and u z = 
0, respectively. The other two solutions diverges as the first and the second power of z, 
respectively, they are referred to as general (non-Bloch) Floquet modes. 

Deviation of the frequency uj from uj removes the triple degeneracy (|75jl of the matrix 
Tl, as seen from Eq. (|7ip. The modified matrix Tl can now be reduced to a diagonal form 
with the set (J5%j) of four eigenvectors comprising two extended and two evanescent Bloch 
solutions. 
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E. Symmetry considerations 

In Section 2, we discussed the relation between the symmetry of the axial dispersion 
relation of a periodic stack, and the phenomenon of AFM. At this point we can prove that 
indeed, the axial spectral asymmetry (J 17)) is a necessary condition for the occurrence of 
the stationary inflection point and for the AFM associated with such a point. As we have 
seen earlier in this Section, the stationary inflection point relates to a triple root of the 
characteristic polynomial F(X) from Eq. ([60)1 . Since F(X) is a polynomial of the fourth 
degree, it cannot have a symmetric pair of triple roots, that would have been the case for 
axially symmetric dispersion relation. Hence, only asymmetric axial dispersion relation u (k) 
can display a stationary inflection point (J13j) or, equivalently, (J 16)1 . as shown in Fig. Efb). In 
this respect, the situation with the AFM is somewhat similar to that of the frozen mode in 
unidirectional magnetic photonic crystals 16]. The difference lies in the physical nature of 
the phenomenon. The bulk spectral asymmetry leading to the effect of electromagnetic 
unidirectionality, essentially requires the presence of nonreciprocal magnetic materials. By 
contrast, the axial spectral asymmetry (|17j) along with the AFM regime can be realized 
in perfectly reciprocal periodic dielectric stacks with symmetric bulk dispersion relation 
()22|) . On the other hand, the axial spectral asymmetry essentially requires an oblique light 
incidence, which is not needed for the bulk spectral asymmetry. 

Another important symmetry consideration is that in the vicinity of the stationary in- 
flection point 1)13)1 ■ all four Bloch eigenmodes (jf)8j) must have the same symmetry, which 
means that all of them must belong to the same one- dimensional irreducible representation 
of the Bloch wave vector group. This condition is certainly met when the direction defined 
by (n x , n y ) is not special in terms of symmetry. Let us see what happens if the above con- 
dition is not in place. Consider the situation 1)31)). when at any given frequency u and fixed 
(n x ,n y ) = (n x ,0), two of the Bloch eigenmodes are TE modes and the other two are TM 
modes. Note, that TE and TM modes belong to different one-dimensional representations 
of the Bloch wave vector group. In such a case, the transfer matrix Tl can be reduced to 
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Tr 



the block-diagonal form 

r T n Tia 

T 2 i T 22 
T 33 T 34 
T 43 T 44 

The respective characteristic polynomial F(X) degenerates into 

F(X) = F TE (X)F TM (X) (83) 

where Fte{X) and Ftm{X) are independent second degree polynomials describing the TE 
and TM spectral branches, respectively. Obviously, in such a situation, the transfer matrix 
cannot have the nontrivial canonical form Eq. (J76j) . and the respective axial dispersion 
relation cannot develop a stationary inflection point (|69J) . regardless of whether or not the 
axial spectral asymmetry is in place. 



IV. THE AFM REGIME IN A SEMI-INFINITE STACK 



A. Boundary conditions 



In vacuum (to the left of semi-infinite slab in Fig. [TJ the electromagnetic field \lV (z) is 
a superposition of the incident and reflected waves 

*v (z) = (z) +V R (z), at z < 

At the slab boundary we have 

m v (0) = *r (0) + * R (0) = $j + $ K . 



where 



E R,y 

Hr, x 



e i,x 
E i,y 

-Ei^n x n y n~ x - E I>y (1 - n 2 x ) n z l 



n y n z 



H 



R,y . 



E I)X (l - n 2 y ) n z 1 + Ei ty n x 

e r,x 
E R,y 



E Rx n x n y n z 1 + E Ry {l-n 2 x )n z l 
-E R , X (l - n 2 ) n~ l - E R ^ y n x n y n- 1 _ 



(84) 



(85) 



(86) 
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The complex vectors Ej, Hf and E R , H R are related to the actual electromagnetic field 
components E/, Hj and E R , H R as 

E/ = e ^{n x x+n y y)^ i ^ ? jj^ = (n x x+n v y) ^ 

E R = e^ {n * x+n y y) E R (z), H R = e^ {n * x+n v y) H R . 

The transmitted wave {z) inside the semi-infinite slab is a superposition of two Bloch 
eigenmodes 

^ T (z) = (z) + ^ 2 (z) , at 00. (87) 

(in the case of a finite slab, all four eigenmodes would contribute to The 
eigenmodes (z) and ^2 {z) can be both extended (with u x > 0), one extended and one 
evanescent (with u x > and ^sk > 0, respectively), or both evanescent (with > 0), 
depending on which of the three cases ()63|). or (JMj) we are dealing with. In particular, 
in the vicinity of the AFM (e.g., the vicinity of ujq in Fig. Efb)), we always have the situation 
(jHSJ). Therefore, in the vicinity of AFM, {z) is a superposition of the extended eigenmode 
\l/ ex . (z) with the group velocity u z > 0, and the evanescent mode ^ ev (z) with > 

m T (z) = m ex (z) + $> ev (z) , at z > 0. (88) 

The asymptotic expressions for the respective wave vectors k ex and k ev in the vicinity of 
AFM are given in Eq. (jUJ) . 

When the frequency u exactly coincides with the frequency o>o of the AFM, the repre- 
sentation (|HHj) for ty T (z) is not valid. In such a case, according to Eq. (|73)l. there is no 
evanescent modes at all. It turns out that at uj = ujq, the electromagnetic field inside the 
slab is a superposition of the extended mode \I/fc ( z ) an d the (non-Bloch) Floquet eigenmode 
^01 {z) from Eq. (JHOJ) 

*r (z) = ^ fco (z) + ^01 (z) , at = cu and 2; > 0. (89) 

Since the extended eigenmode ( z ) has zero axial group velocity u z , it does not contribute 
to the axial energy flux S z . By contrast, the divergent non-Bloch contribution \l/oi {z) is 
associated with the finite axial energy flux S z > 0, although the notion of group velocity 
does not apply here. The detailed analysis is carried out in the next subsection. 

Knowing the eigenmodes inside the slab and using the standard electromagnetic boundary 
conditions 

$ T = $/ + ® R , where $ = ^ (0) , (90) 
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one can express the amplitude and composition of the transmitted wave and reflected 
wave ^r, in terms of the amplitude and polarization of the incident wave This gives 
us the transmittance and reflectance coefficients (j2J) of the semi-infinite slab, as well as the 
electromagnetic field distribution \l/<r (z) inside the slab, as functions of the incident wave 
polarization, the direction n of incidence, and the frequency u. 



B. Field amplitude inside semi-infinite slab 

In what follows we assume that u can be arbitrarily close but not equal to u>o, unless 
otherwise is explicitly stated. This will allow us to treat the transmitted wave {z) as a 
superposition (JHH|) of one extended and one evanescent mode. Since evanescent modes do 
not transfer energy in the z direction, the extended mode is solely responsible for the axial 
energy flux S z 

S z (*r) = S z (*ex) ■ (91) 

According to Eq. (jl()9j) . S z does not depend on z and can be expressed in terms of the 
semi-infinite slab transmittance r from Eq. (J2J) 



S^r^Sjj^rS!, (92) 

where Sj = is the axial energy flux of the incident wave, which is set to be unity. 

The energy density W ex associated with the extended mode ty ex (z) can be expressed in 
terms of the axial component u z of its group velocity and the axial component S z (\l/ ex ) of 
the respective energy density flux 

W ex = u- z l S z = rSj. (93) 

In close proximity of the AFM frequency ujq, we have according to Eq. (fT3*j) 

w-w «k'(fc-^) 3 , (94) 
o 

where u/ " is defined in Eq. ()73|) . Differentiating Eq. (J9~4"j) with respect to k 

%) - \< (k ~ ~ ^ K') 1/3 (u, - -o) 2/3 , (95) 

' kx fey 

and plugging Eq. (j93|) into (|93|) yields 

w ~ ~ w* tSi {uj " yl " {uj ~ w ° r2/s ' (96) 
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where the transmittance r depends on the incident wave polarization, the frequency oj, and 
the direction of incidence (n x ,n y ) = (ck x /uj,cky/uj). Formula implies that the energy 
density W ex and, therefore, the amplitude \^ ex (z)\ = \& ex \ of the extended mode inside the 
stack diverge in the vicinity of the AFM regime 

|$ e ,| ~ \/W^ ~ y^S~i K')~ 1/6 |w - u \- 1/3 as u -> w . (97) 

The divergence of the extended mode amplitude |$ ex | imposes the similar behavior on 
the amplitude |\& ei ,(0)| = \& ev \ of the evanescent mode at the slab boundary Indeed, the 
boundary condition (J9U)) requires that the resulting field $^ — &ex + &ev remains limited to 
match the sum $j + $^ of the incident and reflected waves. The relation (|9U|) together with 
(jHZj) imply that there is a destructive interference of the extended $ ex and evanescent $ e „ 
modes at the stack boundary 

$ ex « -$ e „ « iTv^S; K") _1/6 (w - uJ y 1/3 $ fco as w -> w 0j (98) 

Here $fc is the normalized eigenvector of Tl in Eq. (|7fij) : is a dimensionless parameter. 
The expression (|98j) is in compliance with the earlier made statement (|7Hjl that the column- 
vectors $ ea; and $ et , become collinear as w —>■ ujq- 

1. Space distribution of electromagnetic field in the AFM regime 

The amplitude \^ ex {z)\ of the extended Bloch eigenmode remains constant and equal 
to \<& ex \ from (|97j). while the amplitude of the evanescent contribution to the resulting field 

decays as 

\^ ev {z)\ = \3> ev \ e"* 3 *", where %k ev « ^V 3 (^")- 1/3 |u; - w | 1/3 • (99) 

At 2 ^> (Qkev) -1 , the destructive interference (}9"H)) of the extended and evanescent modes 
becomes ineffective, and the only remaining contribution to \I/x (z) is the extended mode 
ty ex (z) with huge and independent of z amplitude (J9"7)l . This situation is graphically demon- 
strated in Fig. El 

Let us see what happens when the frequency u tends to its critical value ujq. According 
to Eqs. (J5Z!) and ((HHD , at z = NL, N = 0, 1, 2, the resulting field ^ T (z) inside the slab 
can be represented as 

^ T ( z ) = § ex e izkex + ® ev e izke \ (100) 
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Substituting k ex and k ev from Eq. (|73jl in ()100j) . and taking into account the asymptotic 
relation (J98)l . we have 

*' « « (•* + 2/f \S 6 " 3 (s + Tj^j) *") as ^ - ^ (101) 

Although this asymptotic formula is valid only for z = NL, N = 0,1,2,..., it is obvi- 
ously consistent with the expression (J8U|) for the non-Bloch solution \J/ 10 (z) of the Maxwell 
equation (jHHj) at to = loq. 

2. The role of the incident wave polarization 

The incident wave polarization affects the relative contributions of the extended and 
evanescent components to the resulting field ^>t (z) in Eq. (|88j). In addition, it also affects 
the overall transmittance O). The situation here is similar to that of the normal incidence, 

n 

considered in [16||. There are two special cases, merging into a single one as u — > ujq- 
The first one occurs when the elliptic polarization of the incident wave is chosen so that it 
produces a single extended eigenmode ^ ex (z) inside the slab (no evanescent contribution to 
i z ))- ha this case, (z) reduces to \l/ ex (z), and its amplitude \^>t (z)\ remains limited 
and independent of z. As u> approaches u;o, the respective transmittance r vanishes in this 
case, and there is no AFM regime. The second special case is when the elliptic polarization of 
the incident wave is chosen so that it produces a single evanescent eigenmode ^ ev (z) inside 
the slab (no extended contribution to ^t{z))- In such a case, ^t( z ) reduces to "$? ev (z), 
and the amplitude decays exponentially with z in accordance with Eq. The 

respective transmittance r in this latter case is zero regardless of the frequency u, because 
evanescent modes do not transfer energy. Importantly, as u approaches co>o, the polarizations 
of the incident wave that produce either a sole extended or a sole evanescent mode become 
indistinguishable, in accordance with Eq. (J75j) . In the vicinity of the AFM regime, the 
maximal transmittance r is achieved for the incident wave polarization orthogonal to that 
exciting a single extended or evanescent eigenmode inside the semi-infinite stack. 
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C. Tangential energy flux 



So far we have been focusing on the axial electromagnetic field distribution, as well as 
the axial energy flux S z inside semi-infinite slab. At the same time, in and near the AFM 
regime, the overwhelmingly stronger energy flux occurs in the tangential direction. Let us 
take a closer look at this problem. 

The axial energy flux S z is exclusively provided by the extended contribution ^> ex (z) to 
the resulting field ^t(z), because the evanescent mode ^ ev (z) does not contribute to S z . 
Neither \^ ex \ nor S z depends on z (see Eq. (I109jl ). By contrast, both the extended and 
the evanescent modes determine the tangential energy flux S T (z). Besides, according to Eq. 
(jl(J9j) . the tangential energy flux depends on z. Far from the AFM regime, the role of the 
evanescent mode is insignificant, because ^ ev (z) is appreciable only in a narrow region close 
to the slab boundary. But the situation appears quite different near the AFM frequency. 
Indeed, according to Eq. (|9*9*j) . the imaginary part of the respective Bloch wave vector k ev 
becomes infinitesimally small near the critical point. As a consequence, the evanescent mode 
extends deep inside the slab, so does its role in formation of S T (z). The tangential energy 
flux S T (z) as function of z can be directly obtained using formula and the explicit 

expression for (z) = ^ ex (z) + ty ev (z). Although the explicit expression for S T (z) is 
rather complicated and cumbersome, it has very simple and transparent structure. Indeed, 
the tangential energy flux can be represented in the following form 

S T (z) = u T W (z) , 

where the tangential group velocity u T behaves regularly at uj = u . Therefore, the magni- 
tude and the space distribution of the tangential energy flux S T (z) in and near the AFM 
regime literally coincides with that of the electromagnetic energy density W (z), which is 
proportional to ( z )\ 2 - A typical picture of that is shown in Fig. 0a). 

V. SUMMARY 

As we have seen in the previous Section, a distinctive characteristic of the AFM regime 
is that the incident monochromatic radiation turns into a very unusual grazing wave inside 
the slab, as shown schematically in Figs. El and H3 Such a grazing wave is significantly 
different from that occurring in the vicinity of the total internal reflection regime, where the 
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transmitted wave also propagates along the interface. The most obvious difference is that 
near the regime of total internal reflection, the reflectivity approaches unity, which implies 
that the intensity of the transmitted (refracted) wave vanishes. By contrast, in the case of 
AFM the light reflection from the interface can be small, as shown in an example in Fig. 
Ufa). Thus, in the AFM gnificant portion of the incident light gets converted into 

the grazing wave (the AFM) with huge amplitude, compared to that of the incident wave. 
For this reason, the AFM regime can be of great utility in many applications. 

Another distinctive feature of the AFM regime relates to the field distribution inside 
the periodic medium. The electromagnetic field of the AFM can be approximated by a 
divergent Floquet eigenmode \l/i (z) from (|8L)jl. whose magnitude |^ 10 (z)\ 2 increases as z 2 , 
until nonlinear effects or other limiting factors come into play. In fact, the field amplitude 
inside the slab can exceed the amplitude of the incident plane wave by several orders of 
magnitude, depending on the quality of the periodic array, the actual number of the layers, 
and the width of the incident light beam. 

Looking at the z component of light group velocity and energy flux, we see a dramatic 
slowdown of light in the vicinity of the AFM regime, with all possi ble p ractical applications 
extensively discussed in the literature (see, for example, 9|, [lOj, |ll|, |l2], and references 
therein). In principle, there can be a situation when the tangential components (u x , u y ) of the 
group velocity also vanish in the AFM regime, along with the axial component u z . Although 
we did not try to achieve such a situation in our numerical experiments, it is not prohibited 
and might occur if the physical parameters of the periodic array are chosen properly. In 
such a case, the AFM regime reduces to its particular case - the frozen mode regime with 
u = inside the periodic medium. This regime would be similar to that considered in [16| , 
with one important difference: it is not related to the magnetic unidirectionality and, hence, 
there is no need to incorporate nonreciprocal magnetic layers in the periodic array. The 
latter circumstance allows to realize the frozen mode regime at the infrared, optical, and 
even UV frequency range. 
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VI. APPENDIX 1. J-UNITARITY OF THE TRANSFER MATRIX 

Let n x n matrix T (z) satisfies the following Cauchy problem 

d z T(z) =UA(z)T(z) , T(0) = I. (102) 

where A (z) = JM (z) is a Hermitian matrix. Let us prove that the unique solution T (z) 
for Eq. (jlU2|) is a J-unitary operator 

(z) = JT- 1 (z) J' 1 . (103) 
To prove it, notice that Eq. ()102|) implies 

d x T ] (z) = -T f (x) iA (z) J, T f (0) = /. (104) 
Now, let us find the respective Cauchy problem for T^ 1 (z). Since 

d z [T (z) T- 1 (z)] =0 = T(z) [d z T- 1 (z)] + [d z T (z)} T~ l (z) , 

we have 

d.T" 1 (z) = -T- 1 (z) [d z T (z)jr- 1 (z) , 
which in a combination with Eq. (|102j) yields 

d z T~ l (z) = -T- 1 (z) UA (z) , T- 1 (0) = /. (105) 

Finally, multiplying both sides of the equality ()105j) by J and using the fact that J 2 = I, 
we get the following Cauchy problem for JT -1 (z) J 

d z [JT' 1 (z) J]=- [JT' 1 (z) J] iA (z) J, JT' 1 (0) J — I. (106) 

which is identical to that for Tt (z) from Eq. ([104)1 . Since both Cauchy problems (|104|) and 
f!106p have unique solutions, their similarity implies the relation ()1U3|) of J-unitarity. 
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VII. APPENDIX 2. ENERGY DENSITY FLUX 



The real-valued Poynting vector is denned by 

S (f) = — 9? [E* (f) x H (f)] . 
8n 

Substituting the representation (|3*7|) for E (r) and H (r) in Eq. (jl07|) yields 

S(r) = S(z) = E (:) x H (:) 

07T 



(107) 



(108) 



implying that none of the three Cartesian components of the energy density flux S depends 
on the transverse coordinates x and y. Energy conservation argument implies that the 
component S z of the energy flux does not depend on the coordinate z either, while the 
transverse components S x and S y may depend on z. Indeed, in the case of steady-state 
oscillations in a lossless medium we have, with consideration for Eq. 



V • S = d z S z (z) = 
which together with Eq. ()108|) gives 

S z (f) = S z = const, S x (r) = S x (z) , S y (f) = S y (z) . 
The explicit expression for the z component of the energy flux ()108|) is 

s z = \ [e*h v - e;h x + e x h; - e v h:] = ~ {% j*) . 



(109) 



(110) 



The tangential components of the energy flux can also be expressed in terms of the column 
vector \1/ (z) from Eq. ()38|) . Using the expressions (pfUj) for E z and H z and eliminating these 
field components from S (z) in Eq. ()108|) yields 



in) 



where G T and G„ are Hermitian matrices 
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Both G x and G y are functions of the Cartesian coordinate z, frequency u, and the direction 
n of incident wave propagation. 
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